Spin and density overlaps in the frustrated Ising lattice gas 
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We perform large scale simulations of the frustrated Ising lattice gas, a three-dimensional lattice 
model of a structural glass, using the parallel tempering technique. We evaluate the spin and density 
overlap distributions, and the corresponding non-linear susceptibilities, as a function of the chemical 
potential. We then evaluate the relaxation functions of the spin and density self-overlap, and study 
the behavior of the relaxation times. The results suggest that the spin variables undergo a transition 
very similar to the one of the Ising spin glass, while the density variables do not show any sign of 
transition at the same chemical potential. It may be that the density variables undergo a transition 
at a higher chemical potential, inside the phase where the spins are frozen. 



I. INTRODUCTION 

In the last years, a scenario has emerged for the theo- 
retical description of structural glasses, in which the be- 
havior of glasses is characterized by two different temper- 
atures. The higher temperature Tc, identified with the 
critical temperature of the ideal mode-coupling theory 
101 , corresponds to a crossover in the dynamical behavior 
of the glass, characterized by a single relaxation time for 
T > Tc, and by a two-step relaxation for T < Tc- Be- 
low Tc, the free energy landscape in the phase space is 
splitted into an extensive number of valleys, separated by 
high barriers, and the dynamics is separated in a fast mo- 
tion inside the valley, and a much slower motion among 
different valleys. The number of valleys accessible to the 
system at the temperature T, can be expressed by the 
formula JV = exp(7VI](r)), where N is the number of 
particles in the system, and S(T) is the so-called config- 
urational entropy. The second temperature Tk, identi- 
fied with the Kauzmann temperature Q and the critical 
temperature of the Adam and Gibbs theory Q, is char- 
acterized by the vanishing of the configurational entropy 
E(T), and the divergence of the relaxation time of the 
system with a Vogel-Fulcher-Tamman law. 

This picture is corroborated by the analogy with a 
class of mean field models, the p-spin glasses. These 
models undergo a dynamical transition at a tempera- 
ture Td, where the phase space splits into an extensive 
number of metastable states, and the correlation func- 
tions show a singularity of the same type of the one pre- 
dicted by mode-coupling theory Q. At a lower temper- 
ature Ts the number of metastable states becomes non- 
extensive, that is the configurational entropy vanishes, 
and for T < Tg the model shows a 1-step replica symme- 
try breaking. It has been longly debated to what extent 
this analogy, between structural glasses and p-spin mod- 
els, can be pushed forward. In particular one may ask 
if structural glasses exhibit some kind of replica sym- 
metry breaking, and if they are in the same universal- 
ity class of p-spin models or not. There are presently 



some results that point in this direction. First principle 
computations of the equilibrium thermodynamics of sim- 
ple fragile glasses seem to be consistent with this picture 
and the off-equilibrium fluctuation-dissipation ratio, 
in molecular dynamics simulations of supercooled liquids 
shows the same pattern observed in p-spin glasses 
and predicted for 1-step RSB models 10,^. 

Nevertheless, there are still some important differences 
between structural glasses and p-spin models. In p-spin 
glasses the relevant variables are the spin orientations, 
and there is explicit quenched disorder in the Hamilto- 
nian. On the other hand, in glasses the relevant variables 
are the particle positions, and there are no quenched in- 
teractions in the Hamiltonian, but the disorder is "self- 
generated" by the geometrical hindrance that the par- 
ticles exert on each other. Furthermore, p-spin glasses 
are mean-field models, while structural glasses live in fi- 
nite dimension. This means among other things that 
metastable states in p-spin models have infinite lifetime, 
and the dynamical transition Td corresponds to a di- 
vergence of the relaxation times, while in finite dimen- 
sional models the barriers between free energy valleys 
can be overcome by "hopping" processes, which restore 
ergodicity even below the crossover temperature Tc. Fur- 
thermore the p-spin model, if studied in three dimen- 
sions, loses many of its mean field "glassy" properties, 
and shows a transition similar to the full replica symme- 
try breaking transition of Ising spin glasses, though with 
some remarkable differences |9|. 

The frustrated Ising lattice gas was introduced some 
time ago as a simple finite dimensional lattice model of a 
glass- forming liquid pO) , in order to overcome some of the 
limitations of p-spin models. Each lattice site carries two 
kinds of variables, a lattice gas variable rii = 0, 1, which 
represent the presence or absence of a particle on the i- 
th site, and an Ising spin variable S'i, which represent 
an internal degree of freedom of the particle, such as for 
example the orientation of a non-symmetrical molecule. 
The Hamiltonian of the model is 
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H = -/^(l - eijSiSj)ninj ~ fi'^n^, (1) 

where Cij = ±1 are quenched variables. In the hmit 
J — > oo, the first term of the Hamiltonian implies that 
two nearest neighbor sites can be simultaneously occu- 
pied by two particles only if their spin variables satisfy 
the constraint CijSiSj = 1. Therefore, if we identify the 
variables Si with the orientation of a non-symmetrical 
molecule, this condition means that two molecules can 
be near only if their relative orientation is appropriate. 

Being constituted essentially by diffusing particles, this 
model is suited to study quantities like the diffusion coef- 
ficient, or the density autocorrelation functions, that are 
usually important in the study of liquids. Indeed, the 
model has proven to reproduce fairly well many features 
of supercooled glass-forming liquids, as for example the 
"cage effect" . At low temperature and high density, the 
model shows a two step relaxation in the self correlation 
function and in the mean square displacement. Further- 
more, being a finite dimensional model, it is suitable to 
study activated processes, that are absent in mean field 
models. 

Although the model has been extensively studied by 
Monte Carlo simulations ||l^ and in mean field [|ll| , there 
are still many unsolved problems, concerning the type of 
transition presented by the model. In particular, it would 
be interesting to study if, in finite dimension, density 
and spin variables become critical at the same point, and 
what kind of replica symmetry breaking they show in the 
spin glass phase. 

In this paper, we study the static and dynamical equi- 
librium properties of the model in the limit J — > oo. 
We will evaluate the equilibrium overlap distribution of 
spin and density variables, the equilibrium autocorrela- 
tion functions of the self-overlaps, and the self-diffusion 
coefficient. 



II. SPIN AND DENSITY OVERLAP 
DISTRIBUTIONS 

We have simulated the frustrated Ising lattice gas by 
means of the parallel tempering technique . With 

this technique one can thermalize the system at high 
chemical potential (high density) , where the conventional 
Monte Carlo algorithms suffer of extremely long auto- 
correlation times. One has to simulate several identical 
replicas of the system, at different chemical potentials 
/io < • • • < /^n, where /io corresponds to a low density 
and very short autocorrelation time, and /i„ to the high- 
est value of the chemical potential that one wants to 
investigate. Each step of the algorithm consists of the 
following two substeps: for each replica, perform a con- 
ventional Monte Carlo step with the given temperature 
and chemical potential; for each pair of replicas with ad- 
jacent chemical potentials, try to swap them with prob- 
ability Pswap — min(l, cxp(— /3A/iArt)), where An is the 



difference in the number of particles of the two repli- 
cas considered, and A/i the difference of their chemical 
potentials. If one chooses carefully the set of chemical 
potentials, then the replicas will perform a random walk 
over the interval [/io, Mn]- The time needed to go from /iq 
to /i„ and back again is called ergodic time, and can be 
considered as the maximum autocorrelation time of the 
system. 

We have simulated the model for J = oo and /3 = 1, 
between the chemical potentials /imin — 1-69 and /imax — 
10.69, for system sizes 6^, 8^ and 10'^. In the first case 
we have performed the simulation over 12 values of the 
chemical potential, in the second over 16 values, and in 
the third over 25 values. The exact values of the chemical 
potentials were determined by an algorithm which tried 
to optimize the intervals, in order to obtain the same 
swap rate between all the adjacent replicas. For each 
value of the chemical potential we simulated two repli- 
cas, in order to evaluate spin and density overlap. The 
thermalization of the systems was checked by looking at 
the densities of the replicas, and waiting until they did 
not show any sensible drift in time. The thermalization 
time was greater than 10^ steps. 

After having thermalized the systems, at each paral- 
lel tempering step we collected, for each pair of repli- 
cas at the same chemical potential, the spin overlap 
Qs = ji-J2i Sf'^'iSi'^'i and the density overlap qd = 
J! ^f''^i^ where the superscripts a and /? refer to the 
two replicas, and N is the number of spins. We simulated 
the systems for up to 10^ steps for the largest size, and 
checked the symmetry of the spin overlap distribution to 
be sure that the simulation time was sufficient. All the 
relevant quantities were averaged over 32 disorder real- 
izations, and errors were evaluated from the fluctuations 
between different realizations. 

In Fig. [l] the equilibrium distribution of the spin over- 
lap is shown, for the largest size 10'^ and for different 
chemical potentials. For high chemical potential, the 
distribution develops two peaks, separated by a contin- 
uous plateau. This is typical of models with a continu- 
ous replica symmetry breaking, like the Ising spin glass 
in three dimensions. Note that, for the highest chemi- 
cal potential, the distribution between the two peaks is 
not a constant plateau, but is formed by many small 
peaks. This can be due to the fact that, at high chemical 
potential or low temperature, averaging over more than 
32 disorder realizations is needed to obtain a constant 
plateau. The continuous replica symmetry breaking is 
associated with the divergence of the spin glass suscep- 
tibility, defined as XSG — N{Qs)^ where (■ • ■) denotes 
both the thermal average for a given set of interactions 
and the average over the disorder realizations. We have 
evaluated xsG for the different sizes and chemical poten- 
tial, and the result is shown in Fig. ^j. The behavior of 
the susceptibility confirms the presence of a thermody- 
namical second order transition. The exact value of the 
chemical potential at the transition can be evaluated by 
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looking at the Binder parameter g = i(3 — (qi) / {qs)"^) , 
that is shown in Fig. |^. The curves corresponding 
to the different sizes cross at the chemical potential 
/Uc = 3.67. Once we have located the transition, we can 
try to evaluate the critical exponents, using the relation 
Xsg{L,^i) = L'^^^^xsG[L^^''{^^~ fJ-c)], that should be valid 
around the transition, with xsgM a universal curve. In 
Fig. ^ the best fit is shown, obtained for the values = 
and V = 1. Note that these exponents were obtained also 
by Campellone et al. in the three-dimensional version 
of the p-spin model. They are different from the expo- 
nents found in the Ising spin glass, v = 1.7 ± 0.3 and 
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This suggests that the transition 
could belong to a different universality class with respect 
to the Ising spin glass. 

We have then evaluated the equilibrium distribution 
of the density overlap, which is shown in Fig. ^ for the 
largest size and for different chemical potentials. The 
non-linear compressibility — N{{q^) — (qd)'^) is shown 
in Fig. |6| as a function of the chemical potential. The 
arrow in Fig. ^ marks the point where the spin vari- 
ables undergo the transition, as signaled by the crossing 
of the Binder parameter, and the spin glass susceptibil- 
ity diverges. It is evident that this point does not corre- 
spond to a divergence of the non-linear compressibility. 
Therefore the divergence of the spin glass susceptibility 
is due to the fluctuations of spin variables Si, and not 
to the fluctuations of the density variables n^. We have 
also evaluated the Binder parameter of the density over- 
lap (not shown), which shows a non-monotonic behav- 
ior, becoming negative for low density and positive at 
high density, similarly to what is observed in the finite- 
dimensional p-spin 1^. Indeed the curves for different 
sizes do not cross at a definite point, so they cannot be 
used to localize the transition (if any) of the density vari- 
ables. 

It remains to be determined whether or not the density 
variables exhibit a transition not manifested by a diver- 
gence of the non-linear compressibility. Note that for 
very high chemical potential the equilibrium distribution 
of the density overlap develops a secondary minimum. 
This could correspond to a transition of different kind, 
perhaps similar to the 1-step replica symmetry breaking 
transition of the p-spin models. 



III. RELAXATION FUNCTIONS OF SPIN AND 
DENSITY SELF-OVERLAP 

The frustrated Ising lattice gas is known to have 
very large relaxation times at high density or low tem- 
perature [0. Here we want to evaluate the relax- 
ation times of the spin self-overlap, defined as qs{t) — 
jf J2t{Srit')n^(t')Si{t' + t)ni{t' + t)), and the density self- 
overlap qd{t) = ^ EriMt')Mt' + t)) - {nr{t')?), where 
the average (• • •) is done over the reference time t' . We 
have simulated the model for system size 20'^, and 10 



chemical potentials between /x = 2.583 and 3.661, in the 
following manner. We start with an empty system, with 
the interactions randomly chosen, thermalize at the given 
chemical potential for a given time At, save the obtained 
configuration, and then simulate the model saving the 
self-overlaps qs{t) and qd{t) with respect to the config- 
uration at the end of thermalization. Then we repeat 
the process again with a different disorder configuration 
and thermal noise. The thermalization time Ai is at 
least 10 times larger than the time needed to the self- 
overlap to decay to the value 0.1, except for the three 
highest chemical potentials, for which the thermalization 
time was shorter. The self-overlaps were averaged over at 
least 100 different runs, and errors were evaluated from 
the fluctuations between different runs. When evaluating 
the density self-overlap qd{t)^ the quantity J2i{''^ii'^'))^ 
can be taken equal to the average density overlap (qd) 
at equilibrium, as calculated in the parallel tempering 
simulations. 

In Fig. the relaxation functions of the spin self- 
overlap are shown. The solid lines are fit with the 
function ct~^ exp(— (t:/r)^), proposed by Ogielski for the 
Ising spin glass [|l5[. The exponent /3 is nearly constant 
within the errors for all the chemical potentials consid- 
ered, and slightly greater than 0.5, while the exponent 
X varies between 0.2 for the lowest chemical potential to 
0.1 for the highest. The correlation times r are shown in 
Fig. ^, excluding the last three points, which are likely 
to suffer from insufficient thermalization or finite size ef- 
fects. A power law fit r ~ |/i — fJ-d^^'^ gives a dynami- 
cal exponent z = 7.4, slightly greater than the exponent 
z = 6.0 ± 0.8 found by Ogielski for the Ising spin glass. 

In Fig. ^ the relaxation functions of the density self- 
overlap are shown, for the same system size and chemi- 
cal potentials. Note that the relaxation times grow very 
slowly with respect to those of the spin self-overlap. This 
supports the conclusion that the transition at fic = 3.67 
does not involve the density variables rii . The latter prob- 
ably undergo a transition at a higher chemical potential, 
inside the phase where the spin variables are frozen. 



IV. DIFFUSIVITY 

We have simulated the model with a purely diffusive 
dynamics, in the following way. We start with an empty 
lattice, with random interactions, and then slowly raise 
the chemical potential, until a given density is reached. 
Then we switch to a purely diffusive dynamics, with con- 
served number of particles, and thermalize the system at 
the given density. After having thermalized the system, 
we collect the mean square displacement {r'^{t)) of the 
particles as a function of time. The long time regime of 
the mean square displacement is of the form (r^(t)) = Dt, 
from which we extract the diffusion coefficient D. To 
each density, we associate a chemical potential from the 
equilibrium relation between the two quantities. 
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In Fig. |lG the diffusion coefficient D is sliown for a sys- 
tem size 16^^^ as a function of the chemical potential. For 
high chemical potential, it can be well fitted by an Ar- 
rhenius form, D = ae~°'^. Therefore, the diffusion of the 
particles seems to stop only at /i — > oo, which corresponds 
to T ^ 0. The arrow marks the point where the spin vari- 
ables undergo the spin glass-like transition: apparently 
no anomaly in the diffusivity shows up in correspondence 
of the transition. In the inset, the diffusivity as a func- 
tion of the density is shown. Note that for — > oo the 
density goes to a maximum value Pmax — 0.68. The dif- 
fusivity can be well fitted by a power law D — a(po ^ pV i 
where po — 0.681 and 7 = 1.38. 



V. CONCLUSIONS 

We have studied the static and dynamical properties 
of the frustrated Ising lattice gas at equilibrium. A spin 
glass-like transition is found in the spin variables, sig- 
naled by the crossing of the Binder parameter, the diver- 
gence of the non-linear susceptibility, and the develop- 
ment of a continuous replica symmetry breaking in the 
spin overlap distribution. The equilibrium autocorrela- 
tion functions of the spin overlap can be well fitted by an 
Ogielsky form, with a correlation time diverging at the 
critical point. On the other hand, the density variables 
seem to be affected little by the transition, showing no 
divergence either in the non-linear compressibility, or in 
the autocorrelation time. 

The freezing of the model at the chemical potential 
/ic is therefore connected with a second order transition 
in the spin variables, more similar to the freezing of the 
Ising spin glass than to the mode-coupling transition of 
structural glasses. One cannot exclude that the density 
variables undergo a p-spin-like transition at a higher den- 
sity, characterized by a 1-step replica symmetry breaking 
and a discontinuity of the Edwards- Anderson parame- 
ter defined in terms of density variables. This fact is 
suggested by the development of a secondary peak in 
the density overlap distribution at very high chemical 
potential, as well as by the measurements of the off- 
equilibrium fluctuation-dissipation ratio but more 
work is needed to clarify this point. 
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FIG. 2. Spin glass susceptibility xsg as a function of the 
chemical potential, for sizes 6^, 8^ and 10^. 
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FIG. 4. Finite size scaling plot of the spin glass suscepti- 
bility, with /ic = 3.67. The exponents that give the best data 
collapse are r] = and v = 1. 
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FIG. 5. Density overlap distribution P{qa) for size 10^ and FIG. 7. Relaxation functions of the spin self-overlap, for 

chemical potentials /x=3.08, 4.12, 5.79, 10.69. system size 20^ and chemical potentials = 2.583, 2.665, 

2.747, 2.829, 2.911, 2.997, 3.083, 3.264, 3.456, 3.661. Contin- 
uous lines are fits with the function ct~^ exp(— (t/r)''). 
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FIG. 6. Non-linear compressibility function of the 

chemical potential, for sizes 6^, 8^ and 10^. The arrow marks 
the point where spin variables display the transition. 
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FIG. 8. Relaxation times r of the spin self-overlap, as 
obtained by the fits of Fig. ^, for chemical potentials 
2.583 < fJ. < 3.083. The straight line is a fit with the function 
1^ — ^cp^'^, and He — 3.67 fixed, which gives zi/ = 7.4. 
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FIG. 9. Relaxation functions of the density self-overlap, for 
the same system size and chemical potentials of Fig. |^. 
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FIG. 10. Diffusivity as a function of the chemical potential, 
for a system size 16^ . The solid line is a fit with the Arrhenius 
form D = ae""''. Inset: Diffusivity as a function of the den- 
sity. The solid line is a fit with the power law D — a{po — py , 
with po = 0.681 and 7 = 1.38. The arrows mark the point 
where the spin variables display the transition. 
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